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Tokamak 
Ion  beam 
Plasma  heating 


The  recent  availability  of  intense,  space -charge-neutralized,  energetic  ion  beams  haa  led  us  to 
explore  the  possibility  of  using  them  for  heating  tokamak  plasmas.  The  following  key  problems  for 
this  application  are  Investigated:  (I)  beam  Injection  perpendicular  to  J3,  (it)  beam  propagation  across 
^ In  vacuum,  and  (III)  cross-field  beam  propagation  and  energy  deposition  In  plasma.  It  appears 
possible  that,  for  proper  values  of  the  beam  energy,  radius  and  current,  an  Intense  ion  beam  can 
penetrate  to  tbs  interior  of  • tokamak  plasma  and  deposit  its  energy  there. 
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TOKAMAK  PLASMA  HEATING  WITH  INTENSE  PULSED  ION  BEAMS 
I.  INTRCDUCT  ION 

Recently,  sources  of  vary  intense  ion  beans  have  become  available. 
Experimental  advances,1;2  as  well  as  theoretical  understanding9'*  of 
the  production  of  such  beams,  continue  to  proceed  rapidly.  The 
aval lability  of  these  beams  naturally  leads  to  the  proposal  of  uses 
to  which  they  may  be  put.  Several  such  suggestions  have  recently 
been  advanced.7  In  this  paper  we  propose  and  evaluate  the  use  of 
Intense  Ion  beams  for  tokamak  heating. 

Similar  proposals  have  been  made  for  Intense  electron  beams,  but 
the  problem  of  transporting  an  electron  beam  from  an  external  injection 
point  across  magnetic  surfaces  to  the  plasma  Interior  appears  to  be 
difficult  (cf.  Raf.  8 for  progress  on  this  topic). 

Intense  Ion  beams  can  actually  be  charge  neutralized  after 
they  are  produced  by  dragging  along  an  equal  number  of  electrons. 

(In  this  mode  It  Is  possibly  more  accurate  to  talk  of  an  intense 
"plasma  beam”,  but  we  shall  stick  to  the  term  Ion  beam  because  of  the 
way  the  beam  Is  produced.)  As  Is  well-known,8'11  such  a neutralized 
group  of  Iona  and  electrons  can  move  across  a vacuum  magnetic  field, 
essentially  unimpeded.  If  («*plAli)a  » 1 (u>pl  end  are  the  Ion  plasma 
and  gyro  frequencies  In  the  beam),  by  setting  up  surfece  polerlzation 
charges  In  such  a way  that  the  electric  field  In  the  beam  Is  - V x ft. 
This  has  been  verified  experimentally10*11  In  a number  of  studies 
with  plasma  guns  and  laser  produced  plasmas.  Thus  the  propegetlon 

Note:  Manuscript  submitted  September  80, 1976. 


of  an  Ion  baam  from  an  external  beam  production  point  to  tha  adga 
of  a p Laima  may  ba  possible,  although,  aa  va  shall 
saa  later,  there  are  some  important  complicating  effects  which  alter 
this  simple  picture.  The  propagation  of  tha  beam  across  a magnetic 
field  in  a plasma  and  tha  deposition  of  beam  energy  in  the  plasma 
presents  a more  complex  problem,  whose  solution  is  of  great  importance 
for  the  application  proposed  here. 

The  organisation  of  the  rest  of  this  paper  is  as  follows.  In 
Sections  II  and  III  wa  discuss  tha  production  extraction  of  ion 
beau.  In  Section  IV  we  study  beam  propagation  serosa  a magnetic 
field  in  a vacuum.  It  is  shown  that  the  polarization  fields,  in 
addition  to  allowing  beam  transport  across  B,  also  causa  the  beam  to 
expand  along  B • A simple  criterion  for  such  beam  spreading  to  be 
important  is  given.  Section  V discusses  cross-field  propagation 
and  energy  deposition  in  a plasma.  In  Section  VI  ve  consider  the  use 
of  intense  ion  beams  for  plasma  heating  in  tokamaks. 


II.  BEAM  PRODUCTION 


Th«  slop lest  conception  of  an  Ion  diode  accelerator  la  a planar 
anode-cathode  gap  in  which  the  anode  and  cathode  are  capable  of 
emitting  Iona  and  electron*  respectively.  Solution  of  this  problem 
for  space  charge  limited  flow12  shows  that  the  ratio  of  ion  to  electron 
current  density  is  ( j^/j# ) ■ I.85  (Zm^/m^)^.  Since  the  electrons 
and  ions  both  fall  through  the  eame  potential  energy  difference, 
thla  is  also  the  ratio  of  power  delivered  to  the  ion  beam  to  power 
delivered  to  the  electron  beam.  Thua  moat  of  the  energy  goes  to 
electrons,  and  the  ion  energy  efficiency  is  low.  The  basic  problem 
is  thua  to  prevent  energy  from  being  given  to  electrons  so  that  the 
energy  delivered  to  the  gap  goaa  to  producing  the  ion  beam.  Recently, 
two  methoda  of  accoopllahing  this  have  been  demonstrated:  (1)  the 

reflex  trioda,1*3-5  and  (2)  the  magnetically  Insulated  ion  diode.2*3*6 
A schematic  illustration  of  a reflex  triode  is  shown  in  Fig.  1. 

Two  cathodes  are  located  at  x ■ + L on  either  side  of  an  anode  at 
x - 0 . The  anode  material  and  construction  are  such  that  electrons 
easily  pass  through  it.  Assuming  perfect  transmission,  electrons 
would  simply  be  emitted  from  one  cathode  and  collected  at  the  other 
with  aero  velocity.  Thua  no  energy  would  be  given  to  the  electrons. 
Ions  are  extracted  from  both  sides  of  the  anode  surface  and  pass 
through  the  two  cathodes,  as  shown  in  Fig.  1.  In  practice,  the  anode 
is  not  completely  transparent  to  electrons  (e.g.  electron  scattering 
in  the  anode  will  lead  to  a distribution  of  electrons  trapped  in  the 
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gap  potential  wall13),  and  this  allows  Increased  Ion  currents  as 
wall  • email  electron  currant  collected  at  the  anode.  Nevertheless, 
ion  energy  efficiencies14  in  excess  of  60$  have  been  achieved  (for 
both  baasu  )•  Currently  ion  beams  energies  > 1 MV . end  currents  ** 

1 kA/cm2  have  been  produced.  It  is  to  be  expected  that  such  higher 
energies  and  currents  will  be  reported  in  the  near  future.  One  dis- 
advantage of  the  reflex  trlode  is  that  two  beams  propagating  in 
opposite  directions  are  produced.  Thus  one  must  find  a way  of  simul- 
taneous ly  utilizing  both  beams  (or  else  settle  for  one  half  of  the 
possible  energy  efficiency).  Various  schemes  for  injecting  both  beams 
are  possible,  but  will  not  be  discussed  here. 

Figure  2 illustrates  a magnetically  Insulated  ion  diode.  The 
basic  idea  is  that  a magr  *c  field  of  the  proper  strength  is  applied 
parallel  to  the  anode  and  cathode  surfaces  and  prevents  the  flow  of 
electrons  across  the  gap  but  is  not  strong  enough  to  prevent  ion  flow. 
Roughly  speaking,  the  ion  Larmor  radius  (based  on  the  ion  velocity  in 
falling  through  the  applied  gap  voltage)  exceeds  d,  the  gap  spacing, 
but  the  same  quantity  for  electrons  is  less  than  d.  (Of  course  there 
is  soma  modification  of  the  above  statement  due  to  self-consistent 
electron  and  ion  charge  and  currant  density  within  the  gap  [cf.  Raf.  3].) 

Related  to  the  simple  picture  of  magnetically  insulated  ion  diodes 
in  Fig.  2 is  the  recent  discovery19  that  pinched  electron  diodes 
(cf.  Fig.  3)  can  have  high  ion  efficiencies.  Here  the  magnetic 
insulation  is  from  magnetic  fields  that  arise  naturally  rather  than 
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from  externally  applied  magnetic  fields.  The  magnetic  field  created 
by  the  current  In  the  diode  tightly  pinches  the  electrons  to  the  axis 
but  does  not  Influence  the  more  massive  ions.  In  this  way  the  ion 
beam  to  electron  beam  current  ratio  Is  greatly  enhanced,  especially 
for  large  aspect  ratio  diodes  (l.e.  R/d  In  Fig.  3 large).16 

If  one  considers  an  anode-cathode  gap  with  applied  voltage  V mega- 
volts, then,  according  to  the  Langmuir-Child  law,  the  Ion  current 
density  produced  (In  aaperes/cm2)  Is 

JLC  " 60  v®/2/d2  > (D 

(where  d Is  In  centimeters)  If  It  is  assumed  that  no  electrons  are 
present.  In  actual  fact,  the  presence  of  electrons  In  the  gap  allows 
the  possibility  that  the  space  charge  limited  Ion  current  can  be 
enhanced  by  a factor  a(a  2 ^i^LC^*  For  e3Utt,PleJ  reflex  triode 
operation  with  enhancements  as  large  as  ~ 102  has  been  observed.14 


Finally  we  note  a limitation  on  the  current  pulse  time.  Generally, 
plasma  Is  formed  on  the  anode  and  cathode  surfaces.  In  the  absence  of 
a magnetic  field  these  plassus  expand  toward  each  other.  This  results 
in  shorting  out  the  gap  (called  "gap  closure")  at  a time 

tc-d/(Va+vc)  (2) 
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where  v end  v ere  the  enode  end  cethode  plesme  velocities*  For 

e c 

example,  17  velues  of  tc  ~ 40  nenoeec  have  been  measured  for  d « 0.34cm. 
Application  of  e magnetic  field  transverse  to  the  gap  has  been  observed 17 
to  have  the  beneficial  effect  of  greatly  Increasing  tfi  so  that  longer 
current  pulses  would  be  possible  (e.g*  In  magnetically  Insulated  diodes). 


III.  BEAM  INJECTION 
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At  an  example  of  how  an  ion  beam  can  ba  Injected  into  a tokamak, 
consider  the  situation  shown  In  Fig.  4.  A reflex  trlode  is  situated 
in  a guide  tube  which  leads  into  the  side  wall  of  a tokamak.  The 
trlode  and  the  guide  tube  are  subjected  to  a longitudinal  guide  field 
which  curves  to  merge  with  the  main  tokamak  magnetic  field.  The  ion 
beam  emerging  from  the  right  side  of  the  trlode  is  neutralized  by 
electrons  dragged  off  plasma  on  the  outside  surface  of  the  cathode. 
Thus  a neutralized  beam  will  propagate  down  the  drift  tube.  As  the 
beam  reaches  the  end  of  the  guide  tube,  the  field  lines  curve  to 
eventually  be  at  right  anglas  to  the  initial  beam  velocity.  This 
problem  is  similar  to  that  of  plasma  motion  along  a curved  magnetic 
field  analyzed  years  ago  by  Schmidt.9  Schmidt  demonstrated  that  for 
» 02  a polarization  field  will  be  set  up  by  charge  separation 
within  the  beam  due  to  adiabatic  guiding  center  inertial  drifts, 

E - -V  x B 
~o  ^o  K 

where  ia  the  beam  velocity  in  the  guide  tube.  The  quantity  is 
actually  allghtly  leas  than  - ^ 1 S by  an  amount  of  order  O^/uj^  . 

Thus  the  beam  moves  essentially  unimpeded  by  the  magnetic  field. 

In  order  for  the  above  considerations  to  ba  valid  the  beam  density 
must  be  sufficiently  large  that  cu^  » 0*.  From  the  Langmuir -Child 
law  the  beam  density  is 
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a ■ 3 x 1011  a V/d2(cm”3) 


(5) 


where  V la  In  MV  and  d in  cm,  and  we  have  included  the  enhancement 
factor  a introduced  in  Sec.  II.  n can  be  further  increased  by  con- 
verging the  g - field  in  the  guide  tube.  As  will  be  evident  from 
the  example  treated  in  Sec.  VI,  the  condition  » O2  can  be  easily 
achieved . 

One  could  equally  well  use  a pinched  electron  diode15  in  Fig.  4. 
The  applied  longitudinal  field  could  be  small  compared  to  the  self- 
consistent  pinch  field  in  the  diode,  thus  leaving  the  diode  character- 
istics uneffected.  On  the  other  hand,  the  pinch  field  exists  for 
only  a short  time  and  so  does  not  diffuse  through  the  conducting 
electrodes.  Thus,  once  the  ion  beam  emerges  from  the  right  side  of 
the  cathode  it  is  in  the  purely  longitudinal  applied  field  and  is 
not  effected  by  the  pinch  field. 

It  is  less  obvious  how  to  inject  a neutralized  beam  to  travel 
perpendicular  to  the  applied  tokamak  field,  if  a magnetically  insulated 
ion  diode  (Fig.  2)  is  to  be  used.  The  problem  is  that  the  applied 
insulating  field  will  exist  on  both  sides  of  the  cathode  (it  will 
probably  be  in  existence  long  coopered  to  the  g - field  diffusion 
time  through  the  cathode)  and  will  prevent  electrons  from  being  dragged 
off  the  side  of  the  cathode  as  the  ion  beam  emerges.  Thus  the  emergent 
ions  ere  unneutralized  and  will  only  propagate  a distance  of  the  order 
of  their  Larmor  radius,  a distance  too  short  to  be  of  Interest  in  this 
application.  It  may  be  possible  to  get  around  this  limitation,  but 
it  appears  to  be  difficult. 
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IV.  PROPAGATION  OF  THE  BEAM  IN  VACUUM 

As  was  shown  in  the  previous  section,  when  the  Ion  beam  passes 
from  the  guide  pipe  into  the  tokamak  in  which  it  propagates  perpendicular 
to  8,,  s polarization  electric  £ield  is  set  up  which  gives  the  necessary 
S X 6 drift.  Due  to  the  cylindrical  geometry,  outside  the  beam,  this 
polarization  field  has  a component  parallel  to  J$ . This  leads  to  forces 
on  the  beam  such  that  the  beam  tends  to  expand  along  the  magnetic  field. 
As  it  expands,  the  density  decreases  so  propagation  becomes  more 
difficult.  This  expansion  along  the  field  appears  to  be  the  dominant 
effect  which  limits  the  distance  the  beam  can  propagate  in  the  vacuum. 
This  effect  is  demonstrated  below. 

In  the  reference  frame  moving  with  the  initial  beam  velocity, 

V^,  the  beam  acts  like  a dielectric  with  dielectric  constant 

e * e0(l  + assume  that  the  beam  Is  initially  of  circular 

cross-section  with  radius  r ■ a.  Thus  once  the  beam  enters  the  tokamak 

and  is  propagating  (along  y,  say)  perpendicular  to  the  main  magnetic 

field,  J|  “ the  induced  electric  field  can  be  found  by  transferring 

to  a frame  moving  with  V^.  In  this  frame  we  solve  Laplace's  equation 

in  the  plane  perpendicular  to  the  beam  velocity,  V2^  - o,  using 

cylindrical  coordinates (r,  0) where  0 is  the  angle  to  the  x axis.  The 

boundary  conditions  at  the  beam  edge  and  the  condition  at  infinity 

are  continuity  of  i and  cdd/dr  at  r ■ a («  ■ e for  r > a and 

o 

e ■ «Q  Cl  + for  r < a)  and  i -*  - Eq  x at  r -»  ®,  where 

Eo  ■ VqB.  We  find  that  the  induced  z-component  of  electric  field 
outside  the  beam  is 
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Es  - («/t)aEo  sin  20[l  + (ni/a)pl)2]"1  (lw) 

and  the  electric  field  In  the  beam  la  constant  and  purely  In  the  x- 
dlrectlon, 

E “ 5o  + (u,pi^i)2]'  (**> 


The  polarization  surface  charge  density  at  r - a can  be  calculated 
from  the  difference  between  the  normal  fields  Inside  and  outside  the 
beam.  The  result  Is 


Pg(e)  -2«o 


_eL 


CD2 

pl- 


an2 

d 


E cos 
o 


(5) 


From  the  surface  charge  density,  the  force  per  unit  area  in  the  z 

direction  can  be  calculated  from  F (9)  ■ E (r  ■ a+)  p (0)  to  be 

z z s 


Mote  that  Fz  Is  directed  outward  from  the  beam,  Fz(0)  > 0 for  n > 6 > 0 
and  F (0 ) <0  for  2 tt  > 0 > n. 

Z 

The  way  In  which  the  beam  will  expand  along  the  field  line  due 
to  the  presence  of  this  force  Is  not  at  all  obvious.  Probably  the 
electron  polarization  charge  layer  will  initially  expand  ouch 


faster  than  the  Ions  in  the  positive  charge  layer.  The  electrons  will 
only  be  arrested  when  sufficient  positive  charge  is  left  behind 
to  reverse  the  sign  of  the  parallel  force.  The  electron  sheath  may 
expand  along  a fiald  line  for  a distance  comparable  to  the  radius  of 
the  beam  so  that  the  charge  layer  may  become  grossly  distorted. 

As  a very  crude  approximation,  we  will  say  that  the  beam  expands 
outward  retaining  an  elliptical  shape  of  major  radius  b and  minor 
radius  a equal  to  the  Initial  radius . Then,  integrating  Eq . (6 ) 
(assuming  uj^  » ) over  6 to  determine  the  total  outward  force, 

we  find  an  approximate  equation  for  b 


d2 

"dP" 


~1T  noM 


k 


a«  E2 
o o 


(7) 


where  nQ  is  the  initial  beam  density.  From  (7)  a rough  expression 
for  the  time  for  the  beam  to  expand  along  z by  its  own  radius  is  then 
given  by 


While  the  beam  propagates  forward  at  spaed  V,  the  positive  and 
negative  sheaths  on  the  top  and  bottom  only  propagate  forward  at 
about  V/2  since  the  electric  field  is  reduced  to  approximately  zero 
at  the  outer  edge  of  the  sheaths.8  Thus  the  sheaths  are  constantly 
falling  behind  the  beam  and  are  being  replenished  from  the  interior. 
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From  Che  expression  for  Che  surface  charge  density,  we  find  that 
the  thickness  of  the  sheath  Is  given  by 

2c 

Ls  - VB  cos  0 . (9) 

If  the  sheath  falls  behind  at  a velocity  V/2,  one  can  calculate  that 
the  tine  for  the  beam  to  lose  half  Its  particles  due  to  sheath  loss 
Is  given  roughly  by 

T.  ■ v-  V Tb  (10) 

Cl 

where  t is  the  tine  duration  of  the  bean  pulse.  Clearly,  since 

D 

tupi  » i for  cases  of  practical  Interest,  we  see  that  T.  » Te 
so  that  expansion  along  the  field  lines  Is  the  dominant  effect  limiting 
the  propagation  of  the  beam  In  vacuum. 

Another  possible  physical  effect  which  can  stop  the  beam  propagation 
In  the  vacuum  is  short  circuiting  along  a field  line.  That  is.  If  the 
sane  magnetic  field  line  passes  through  both  the  positive  and  negative 
polarization  charge  layer,  current  can  flow  along  the  field  line  and 
short  circuit  the  charge  separation.  However,  as  was  pointed  out  In 
similar  work  on  octupoles,10  the  distance  along  the  field  line  between 
the  positive  and  negatlva  charge  layers  must  be  short  (not  such  longer 
than  tha  radius  of  the  beam)  In  order  for  this  to  happen.  This  is 
quite  reasonable  Intuitively.  If  an  electron,  after  leaving  the  charge 


12 


I 


I 

layer,  has  traveled  a long  distance  and  has  not  reached  the  positive 
layer,  it  leaves  behind  positive  charge  which  will  pull  it  back  along 
the  field  line. 

For  excels,  in  a tokamak  the  magnetic  field  lines  fill  nested 
magnetic  surfaces  so  that  an  electron  traveling  along  a field  line 
which  passes  through  the  negative  charge  layer  will  ultimately  pass 
through  the  positive  charge  layer;  and  hence  it  might  be  thought  that 
the  field  could  be  shorted  out.  This  consideration,  however,  is 
Irrelevant  not  only  because  the  distance  an  elactron  must  travel  is 
much  longer  than  a beam  radius,  but  also  because  2ttr/c  is  typically 
longer  than  the  duration  of  the  beam  I 


V.  PROPAGATION  AND  ENERGY  DEPOSITION  IN  THE  PLASMA 

After  the  beam  has  propsgstsd  from  ths  wall  through  th«  vacuum 

ragioa  of  ths  tokamak,  It  enters  the  plasma.  We  now  discuss  qualitatively 

how  the  beam  can  propagate  and  deposit  its  energy.  To  begin,  assume  a 

two  dimensional  picture  with  no  variation  along  B.  Since  u>2,  » ft?, 

^ pl  1 

the  electric  field  in  the  beam  must  be  small  in  the  frame  moving  with 
the  initial  beam  velocity,  y^,  (the  beam  acts  like  a perfect  conductor) 
so  that  the  surface  of  the  beam  is  approximate ly  an  aquipotentlal . Since 
the  plasma,  when  it  Jg  x ft  drifts  remains  along  equlpotentlals,  the  plasma 
does  not  interpenetrate  the  beam  but  rather  is  pushed  out  of  its  way. 

Such  a two  dimensional  flow  pattern  in  the  frame  moving  with  the  beam 
is  shown  in  Fig.  5 in  the  plane  perpendicular  to  B.  Of  course  the 
geometry  is  not  two  dimensional,  but  cylindrical.  Each  magnetic  field 
line  in  the  beam  must  be  at  a different  potential  in  order  to  maintain 
the  polarization  electric  field  which  convects  the  beam.  However  the 
potential  of  the  background  plasma  is  the  same  on  different  field  lines. 
Therefore  the  background  plasma  attempts  to  short  circuit  the  polarization 
field  and  stop  the  beam  at  the  outside  of  the  tokamak  plasma,  rather  than 
the  center.  At  the  present  time  the  dynamics  of  the  short  circuiting 
process  is  not  clear  to  us.  However,  even  if  short  circuiting  occurs 
a sufficiently  Intense  beam  will  propagate  through  the  plasma.  To  see 
this  imagine  that  the  polarization  field  is  completely  shorted  out  so 
that  the  forward  motion  of  the  beam  is  converted  to  cyclotron  motion.  Then 
the  diamagnetic  current  on  the  beam  surface  can  exclude  the  tokamak 
magnetic  field  as  long  as  the  tokamak  magnetic  field  energy  density  is 
smaller  than  the  beam  energy  density,  or  that  the  (unneutralized)  beam 


currant  la  large  enough  to  give  field  exclusion.  One  can  easily  show 
that  this  reduces  to 


V /V.  > 1,  (11) 

O A 

where  Is  the  Alfven  speed  in  the  beam.  If  Eq.  (11)  Is  satisfied, 
the  beam  will  exclude  tho  magnetic  field  as  shown  in  Fig.  6 and 


penetrate  the  plasma.  In  order  to  show  the  existence  of  such  a mode 
of  beam  propagation  with  magnetic  field  exclusion  in  the  beam  Interior, 


the  Appendix  demonstrates  a self-consistent  beam  equilibrium  with  the 
above  properties. 

Such  a flow  pattern  shown  In  Figs.  5 and  6 Is  unstable  to  the 
Ka lvln-He Imho ltz  Instability  as  long  as  the  beam  velocity  is  less  than 
a/2  where  Is  the  Alfven  speed.18  [This  is  analogous  to  the  result 
for  the  Ion- Ion  Instability18  where  the  plasma  and  beam  do  interpenetrate. 
There  the  beam  velocity  must  be  less  than  twice  the  Alfven  speed  for 
Instability  to  be  possible.]  As  long  a V <2VA,the  special  growth 
length  of  the  Kelvln-Helmholtz  Instability  Is  given  by18 


k 


(vS 

\ + % 


where  and  pp  are  the  mass  densities  of  the  beam  and  plasma  respec- 
tively. Assuming  that  the  wave  number  which  gives  most  effective 
mixing  is  k ■ a"1,  we  find  1^  **  a(Pb  + p^J/fp^P^)  . Assuming  that 
about  five  growth  lengths  are  needed  to  mix  the  beam  and  plasma,  one 


can  choose  the  beam  radius  and  density  so  that 
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9 


(13) 


•V  ^ pb  PP  ^ r $ A V^p  ( r )3  - 5 

o 

*^*r®  rQ  **dius  of  the  pie  erne . This  la  an  approximate  condition 

for  the  beam  to  propagate  to  the  plasma  interior  and  deposit  its 
energy  there.  If  the  beam  and  central  plasma  densities  are  equal, 
then  from  (13)  the  beam  radius  should  be  about  one  tenth  of  the  plasma 
radius,  where  we  assume  a plasma  density  profile  of  the  form 

pp(r)  " pp(°)  Cl  - (*/*„)*]. 
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VI.  TOKAMAK  HEATING:  AN  EXAMPLE 

We  now  discuss  tbs  problem  of  hosting  s rose  tor 20  to  s tespereture 
of  10  kov  with  lntsnso  ion  bosom.  Assuming  s pulse  time  r,  the  totsl 
energy  of  the  boom  pulse  is 

Eb  - VIT  - 2 * ID-10  tt eaV3/2Tno 

in  Megsjoules.  If  the  tokamak  has  aejor  radius  R end  the  plssas  has 
rsdius  rp  (both  in  cm),  end  average  density  a (in  an*3),  then  the  energy  in 
Megsjoules  to  heat  the  electrons  snd  ions  to  s teqpersture  of  10  kev 

is 

c - 6.4  x 10-21  n tt2  r^. 

P P 

Taking  a voltage  of  5 Megavolt,  (rp/a)  - 10,  n/nQ  - 0.5  and  t - J x 10"8 

sec,  we  find  that  about  ten  pulses  are  needed  for  a reactor  with  a 

radius  of  5 meters.  If  the  central  plasma  (and  beam)  density  is  5 x 1013cm~3, 

and  the  magnetic  field  is  5 x 10 4 Gauss,  then  Eqs.  (11)  and  (13)  are 

satisfied,  and  V < a/2  . Thus  the  beam  should  deposit  its  energy 

near  the  center  of  the  plasms.  Using  Eq.  (3)  we  obtain  for  this 

~ i 

example  d < 0.1  er  cm.  Finally,  let  us  note  that  for  ion  beam  heating 
of  tokamaks  there  ere  many  possible  parameters  to  vary  including  beam 
voltage,  current,  pulse  time,  beam  radius  and  injection  angle.  This 
flexibility  increases  the  likelihood  that  parameters  can  be  found  for 
which  the  beam  will  propagate  to  the  center  of  the  plasma  and  deposit 
its  energy  there. 
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APPENDIX:  Equilibria  for  Croaa  Field  Beam  Propegation  in  Pla8ma 

In  Chla  appendix  we  will  show  that  it  la  possible  to  produce  a 
magnetically  shielded  region  through  which  the  beam  can  propagate. 

Our  interpretation  is  that  the  head  of  the  beam  meets  the  plasma  at 
one  point  and  stops.  In  stopping,  it  sets  up  a diamagnetic  current 
which  cancels  the  ambient  tokamak  field,  allowing  the  rest  of  the  beam 
to  propagate  through.  Thus  the  beam  front  is  continuously  peeled  off, 
creating  a field  free  region  behind  it  through  which  the  remainder  of 
the  beam  propagates. 

For  simplicity,  consider  a two  dimensional  situation  in  the  (x,y) 
plane  perpendiculer  to  Bfy)^,  with  the  beam  propagating  in  x.  Since 
background  plasma  is  assumed  to  short  out  any  electric  fields,  we  set 
E * 0 in  the  laboratory  frame.  The  beam  ions  are  assumed  to  satisfy 
the  Vlasov  equation  so  that  the  beam  ion  distribution  function  depends 
only  on  constants  of  the  motion,  F.  - F. (v  2,  v , P ) where 
Vj2  - vx2  + vy2,  Px  - m^  - eA(y),  and  B(y)  * dA/dy.  Since  the 
electrons  are  much  less  massive  we  use  the  fluid  equations.  From  the 
y - component  of  the  beam  electron  momentum  equation  we  have  -dp^/dy 
- J#  dA/dy  ■ 0 or  ■ -dp^/dA,  where  is  the  beam  electron  current 
(which  la  in  the  x direction)  and  p#  is  the  electron  pressure.  Now 
using  the  V x B Maxwell  equation  we  have 

" ^O  [J1(A)  ■ -dT“]  ' (A1) 
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where 

Jt(A)  - e / vx  d3v, 

which  cen  be  written  J^A)  - - dp^  /dA,  after  ualng  a parts  Integration, 
where  the  yy  beam  ion  pressure  is 

Pt  (A)  - F±  d3v. 

Now  multiplying  Eq.  (Al)  by  dA/dy  and  integrating  over  y we  obtain 

2 57  (^)  + Pe(A)  +Pi(A)  " “2S7Bo2  ' W 

where  p and  p,  go  to  zero  away  from  the  beam,  B is  the  external 
magnetic  field  outside  the  beam.  Equation  (A2)  is,  of  course,  just 
the  y-component  of  the  total  pressure  equation. 

We  now  demonstrate  how  a magnetically  shielded  equilibrium  can 
be  formed.  It  seems  reasonable  to  assume  that  the  beam  electron  « 

random  velocity  is  roughly  the  same  as  the  beam  ion  random  velocity, 
thus  the  electron  pressure  is  much  less  than  the  ion  pressure.  We 
will  begin  by  neglecting  the  electron  pressure.  However,  as  we  will 
shortly  see,  a small  electron  pressure  is  necessary  for  the  formation 
of  the  magnetically  shielded  equilibrium. 

For  concreteness,  let  us  pick  an  ion  distribution  function: 

Ft  - No6 (vA-vxe )6(vx  - §A)6(vf).  ( A3 ) 
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Then  If  electron  pressure  Is  neglected,  Eq.  (A2)  reduces  to 


I (£-)'  ♦ (‘  - (-^tT)4  ■ 5 


(A4) 


Equation  (A4)  Is  analogous  to  the  equation  of  a particle  in  a potential 
'P(A)  • M^NMv^  Q fi  "(  l shown  *n  8oll-d  curve  of  Fig.  7>  where 


A is  analogous  to  the  particle  position  and  y to  time.  There  are  two 

B2 

types  of  solution:  class  1 where  ^ — > N Mvf  , and  class  2 where 

d\L  o ±o 
o 

g2 

~2u~  < Nq^o  * T*ie  vector  potential  and  magnetic  field  are  shown 
as  a function  of  y in  Fig.  8.  Notice  that  class  1 solutions  are  the 
only  reasonable  solutions  for  magnetically  shielded  beam  propagation 
since  the  magnetic  field  does  not  change  sign  on  crossing  the  beam. 

Also  since  B(+  «)  - B(-  ®),  there  is  no  net  ion  current  in  the  x direction. 
However,  class  1 solutions  do  not  completely  exclude  the  field  unless 
B2  p 

-* — “ N Mvf  , in  which  case  the  zero  of  B is  infinite lv  far  from  the 
2|i  o J-o 
o 

position  where  B * Bq.  Notice  that  the  condition  for  field  exclusion 

is  simply  that  v,  > V.  as  stated  in  Sec.  V. 

AO  A 

In  order  for  B(y)  to  go  to  zero  in  a finite  space  the  potential 
«p(A)  must  have  a sharp,  rather  than  smooth,  top  (discontinuous  dA/dy) 
as  shown  by  the  dashed  curve  in  Fig.  J.  Since  the  electron  pressure 
is  very  small  compared  to  the  ion  pressure,  and  the  electron  Larmor 
radius  is  roughly  m/M  compared  to  the  ion  Larmor  radius,  then  the 
electron  pressure  p (A)  can  reasonably  be  considered  to  have  a sharp 


maximum  as  in  Fig.  J 


To  go  back  to  the  description  of  A(y)  via  the  analogy  with  particle 
motion  in  a potential  field,  if  the  particle  energy  is  equal  to  the 
maximum  potential  energy,  then  the  particle  comes  to  rest  on  top  in  a 
finite  time.  It  can  then  remain  at  rest  for  any  arbitrary  time,  and 
finally  roll  off  the  other  side.  The  vector  potential  and  magnetic 
field  corresponding  to  such  a solution  of  Eq.  (A2)  is  shown  in  Fig. 

9 • Notice  that  there  is  an  arbitrary  sized  region  having  B ■ 0 
through  which  the  remainder  of  the  beam  can  propagate. 

Finally,  we  should  point  out  that  the  nature  of  the  equilibrium 
we  have  calculated  is  independent  of  the  particular  form  of  F^.  As 
long  as  p. (A)  has  a maximum  point  somewhere,  and  p (A)  has  a sharp 
top  at  the  maximum  of  p^(A),  such  a magnetically  shielded  equilibrium 
can  form. 
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ANODE  CATHODE 

(x=OI  (x  = d) 

Fig.  2 — Magnetically  insulated  ion  diode.  Electrons  are  confmed  by  the  magnetic 
field  to  the  cross-hatched  sheath  region  adjacent  to  the  cathode 
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ANODE  CATHODE 

Fig.  3 — Electron  pinch  diode,  (bowing  typical  electron  orbits  in  the  self-consistent 
azimuthal  magnetic  field  and  typical  ion  orbits  (which  are  essentially  straight) 


Fig.  4 — Magnetic  field  lines  in  guide  tube  and  near  tokamak  wall. 
A reflex  triode  is  shown  schematically  in  the  guide  tube. 
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Fig.  7 — <fi(A)  versus  A.  Solid  line  is  contribution  of  ions  alone. 
Dashed  line  illustrates  the  effect  of  electrons. 
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CLASS  I SOLUTION 


CLASS  2 SOLUTION 


Fig.  8 — Two  classes  of  solution  for  Eq.  (A4)  (with  ion  pressure  only) 
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